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ABSTRACT 



> 

'• An extremely high-resolution (> 10^) high-S/N (> 10^) solar spectrum has 

been used to measure 15 very weak first overtone {Av = 2) infrared OH lines, 



O ! resulting in a low solar abundance of Aq ~ 8.6 when MARCS, 3D, and spatially 

Q ! and temporally averaged 3D model atmospheres are used. A higher abundance 

r-| ! is obtained with Kurucz [Aq ~ 8.7) and Holweger & Miiller [Aq ~ 8.8) model 

^! atmospheres. The low solar oxygen abundance obtained in this work is in good 

O I agreement with a recent 3D analysis of [O I], O I, OH fundamental {Av = 1) 

^ ', vibration- rotation and OH pure rotation lines (Asplund et al. 2004). The present 

result brings further support for a low solar metallicity, and although using a low 
solar abundance with OPAL opacities ruins the agreement between the calculated 
^ I and the helioseismic measurement of the depth of the solar convection zone, recent 

5^ I results from the OP project show that the opacities near the base of the solar 

convection zone are larger than previously thought, bringing further confidence 
for a low solar oxygen abundance. 

Subject headings: Sun: abundances - Sun: photosphere 



^Affiliated with the Seminario Permanente de Astronomia y Ciencias Espaciales of the Universidad Na- 
cional Mayor de San Marcos, Peru. 
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1. Introduction 

For many years, the accepted solar oxygen abundance was Aq^ ~ 8.9. This value was 
supported by the impressive agreement between the abundance obtained from the [O I] lines 
(Lambert 1978, Aq = 8.92 dex), and the pure rotation (Sauval et al. 1984, Aq = 8.91 dex) 
and fundamental vibration- rotation OH lines (Grevesse et al. 1984, Aq — 8.93 dex). That 
solar oxygen abundance {Aq — 8.9) is high when compared to the abundance of nearby 
young OB stars (e.g. Daflon et al. 2003, Aq — 8.6-8.7 dex). Considering that the Sun is 
about 4.6 Gyr old, its oxygen abundance is supposed to be lower or equal than the present 
day interstellar medium (ISM). The "high" {Aq — 8.9) solar oxygen abundance with respect 
to the ISM could be explained by recent infall of metal-poor material into the local ISM 
(Meyer et al. 1994) or by migration (Wielen et al. 1996) of the Sun from an inner (more 
metal-rich) birthplace in the Galaxy to its present position. 

Another alternative is to invoke errors in the oxygen abundance derived in the Sun (or 
in OB stars). The solar oxygen abundance has been recently revised (Allende Prieto et al. 
2001; Asplund et al. 2004; Allende Prieto et al. 2004; hereafter APOl, A04 and AP04, 
respectively), yielding an abundance about 0.2-0.3 dex lower than previous determinations 
(Sauval et al. 1984; Grevesse et al. 1984; Grevesse & Sauval 1998). Although this result 
solves the apparent oxygen overabundance of the Sun with respect to the local ISM, it poses 
a new problem. The computed base of the solar convective zone (using the low value of 
the solar oxygen abundance) is then shallower than the value obtained from helioseismology 
(BahcaU & Pinsonneault 2004; Bahcall et al. 2004a; Basu & Antia 2004, Turck-Chize et al. 
2004), This new solar problem has been called by Bahcall et al. (2004b) "the convective zone 
problem" . Furthermore, when low solar abundances are employed, the computed sound speed 
radial profile disagrees with sound speeds determined from helioseismological observations 
(Bahcall et al. 2004a; Basu & Antia 2004, Turck-Chize et al. 2004). Therefore, it is 
important to verify the proposed low solar oxygen abundance with other spectral features. 

In the most comprehensive study utilizing difi^erent set of lines available for oxygen 
abundance determination in the Sun, namely the [O I] (0.6 /xm), O I (0.6 - 0.9 /xm), funda- 
mental {Av — 1) vibration-rotation OH (3-4 /xm), and pure rotation OH (9-13 fim) hues, A04 
derived Aq — 8.66 ± 0.05 dex by using a 3D hydrodynamical model of the solar atmosphere. 
The first overtone {Av — 2) OH lines (1.5-1.8 /xm) were probably not included in that work 
because they are so weak that it is very difficult to identify and measure these lines in the 
solar spectrum^. Note that in a recent 3D NLTE analysis of center-to-hmb observations 



^Ax = log {Nx/Nh) + 12 

^After receiving the first report from the referee, we noted that first overtone OH Hnes in the solar 
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of the O I infrared triplet (0.777 //m), AP04 have obtained Aq = 8.70 ± 0.04 dex, this is 
somewhat higher than the abundance obtained from permitted hnes by A04 [Aq = 8.64). 
Considering the uncertainties in the NLTE modehng, both results are equally plausible. 

In this work we present measurements of extremely weak first overtone infrared OH 
lines, which we use for a new determination of the oxygen abundance in the Sun. 

2. Observational Data 

The solar spectrum we use was obtained with the 1-m FTS at the McMarth- Pierce 
telescope (Kitt Peak) by Wallace & Livingston (2003, hereafter WL03). The center-disk 
solar spectrum was obtained at very high spectral resolution (R = 700 000), and corrected 
for telluric hnes (WL03). A portion of this spectrum is shown in Fig. 1. 

The S/N of the spectrum is very high, it typically ranges from 900 to 1700, with a mean 
value of S/N = 1350, as measured from 20 continuum regions in the H band (1.5 - 1.8 /xm). 
This extremely low-noise spectrum allows the measurement of features as weak as a few 0.1% 
of the continuum. 

The printed version of the solar atlas of WL03 provides identification for more than 
seventy OH hnes of the Av — 2 sequence, from the bands (2,0), (3,1), and (4,2). The 
solar atlas of WL03 provides also identification for some features that are blended with the 
infrared OH lines. For the less obvious blends, a line list of atomic and CN lines (Melendez 
& Barbuy 1999) was employed to avoid OH lines that are severely blended. Also, a search 
for new OH features present in the solar spectrum was undertaken using the list of OH lines 
by Goldman et al. (1998, hereafter G98). After careful examination, only 13 of the OH 
lines identified by WL03 were selected as suitable for an abundance analysis; the other lines 
were discarded because of blending with known atomic or CN lines (blends with unidentified 
absorption features were also taken into account, considering only OH lines with FWHM^ = 
0.245 ± 0.025 A). Two additional infrared OH hnes were identified in the process. 

The equivalent widths of the OH lines are shown in Table 1, they were measured by 
fitting gaussian profiles with IRAF. The error in the measurement of these weak features 
is about 10-20 % (0.04-0.08 dex), as estimated from several trials using different continua. 



spectrum were first analyzed by Bonnell & Bell (1993), in a paper on "The gravities of K giant stars 
determined from [O I] and OH features" 

^Bonnell & Bell (1993) reported an average FWHM = 0.093 cm"! for 18 OH solar lines of the Av = 2 
sequence. This is equivalent to FWHM = 0.249 A , in excellent agreement with our results (0.245 A) 
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These errors are consistent with an estimate obtained by using Cayrel's (1988) formula (7), 
that results in 5Wx = 0.11 mA , which for a typical 1-mA OH line corresponds to an error 
of 0.045 dex. 

The errors in Wx could be also estimated from a comparison with measurements by 
Bonnell & Bell (1993), who used the infrared solar spectrum of Delbouille & Roland (1981), 
that has a lower spectral resolution by a factor of 2.33, but with higher S/N by a factor 
of 2.96, therefore both spectra have similar figure of merit F = (R x [S/N])/A (Norris ct 
al. 2001), although note that the Delbouille & Roland solar atlas is not corrected of telluric 
blends. Bonnell & Bell (1993) give Wx for 18 first-ovcrtonc OH lines in the H band, and eight 
of them are in common with our work. They estimated an uncertainty of 9% (0.04 dex) for 
their measurements, but the standard deviation of the oxygen abundance obtained by them 
is 0.09 dex (23%), therefore the errors of their Wx are probably similar to our measurements 
(10-20%). A comparison for the eight OH lines in common is shown in Fig. 2. The mean 
difference (this work - theirs) is 0.07 mA (a — 0.19 mA), equivalent to a mean difference of 
8% ((7= 21%), or 0.03 dex {a = 0.08 dex). 



3. Analysis 

The hst of OH lines and their molecular parameters has been previously described in 
Melendez et al. (2001). It has been employed to determine the oxygen abundance in halo 
(Melendez ct al. 2001) and bulge stars (Melendez et al. 2003). The (^/-values of the first 
overtone OH vibration-rotation lines (X^H) were computed from the Einstein A coefficients 
given by G98, that are the most accurate available in the literature. Moreover, the transition 
probabilities of G98 are in good agreement with previous research (e.g. Holtzclaw et al. 
1993), although significant differences exist with respect to the earliest works (e.g. Mies 
1974). G98 and Holtzclaw et al. (1993) determined Einstein A coefficients employing the 
experimentally derived dipole moment function (DMF) by Nelson et al. (1990, hereafter 
N90). N90 also determined A coefficients, but only for OH lines with low J (J" < 14.5) 
numbers (calculations for somewhat higher J numbers were made available by D. D. Nelson, 
Jr. 1999, private communication). Holtzclaw et al. (1993) noted that for lines of high J 
(and high v), an extrapolation of Nelson's DMF was necessary, although for the lines used 
in the present work the differences due to this extrapolation are negligible. In a different 
approach, G98 used a theoretical DMF outside the range of validity of Nelson's DMF, this 
allowed G98 to extend the calculation of Einstein coefficients to lines with very high J and 
V. Note, however, that for the typical first overtone OH lines observed in the solar spectrum 
(low V and J), any of these approaches give essentially the same result. G98 provides a 
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comparison with the values obtained by Holtzclaw and collaborators, and the agreement is 
excellent. In Fig. 3c it is shown the difference between the log gf value obtained by G98 
and N90; the error bars were derived by N90, estimating the 95% confidence limit in its 
experimental DMF, taking into account both random and potential systematic errors. For 
the lines employed here, the error in the 5f/-values is about 0.03 dex. 

Line positions and excitation potentials are from G98, they were calculated using the 
best available laboratory molecular constants, and the line positions were checked with lab- 
oratory wavenumbers (G98). The adopted value for the dissociation potential is 4.392 eV 
(Huber & Herzberg 1979). The parameters (wavelength, excitation potential, gf -value, 
quantum numbers) of the set of lines used in this work are giving in Table 1. 

The oxygen abundance was obtained using the 2002 version of the program MOOG 
(Sneden 1973) and four model atmospheres: (i) a Kurucz model^ including convective over- 
shooting, although a model with no overshooting (NOVER, Castelli & Kurucz 2004) was 
also used, but the difference in the oxygen abundance is only 0.01 dex (lower with the 
NOVER model); (ii) the temporally and spatially averaged 3D solar model ^ (hereafter 
<3D> model) by A04; (iii) a MARCS solar model (Asplund et al. 1997), (iv) the Holweger 
& Miiller (1974, hereafter HM74) model atmosphere. Additionally, M. Asplund (private 
communication, 2004) have kindly made 3D calculations for a subset of the OH lines. 

A microturbulence of 1 km s~^ was adopted for the ID modeling, yet the lines are so 
weak that the results do not depend on the adopted microturbulence. The calculations were 
done in LTE; this seem to be a good approximation for vibration-rotation OH hues (Hinkle 
& Lambert 1975). 

The resulting abundances obtained with the <3D> model are shown in Fig. 4 as a 
function of excitation potential, reduced equivalent width (Wx/X) and wavelength. As can 
be seen, the derived oxygen abundances show no trends with either excitation potential or 
wavelength. There is an apparent slight trend with reduced equivalent width, but within the 
measurement errors (§2) this is not significant. The use of other model atmospheres result 
in similar plots. 

The mean value obtained from the 15 OH lines and the <3D> model is Aq = 8.59 
dex (cr = 0.06 dex). Full 3D calculations were kindly performed by M. Asplund for 5 lines 
of the A V = 2 sequence with Xexc = 0.30 - 1.24 eV. After rescaling his results for the 
15 OH lines, the full 3D result is 0.02 dex lower than the abundance obtained with the 



^http: / /kurucz. harvard.edu 

^The use of this model was suggested by the referee 
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<3D> model. The MARCS, Kurucz and HM74 models, give Aq = 8.61, 8.71 and 8.80 dex, 
respectively. Our results are in good agreement (except the result with the HM74 model) 
with 3D hydrodynamical simulations of other spectral features by APOl, A04 and AP04, 
who derived Aq = 8.69 ± 0.05 dex, 8.66 ± 0.05 dex, and 8.70 ± 0.04 dex, respectively. 

The main uncertainty in our results is the statistical error, due to the weakness of the 
OH lines. The only previous abundance analysis of the A v = 2 sequence was performed 
by Bonnell & Bell (1993). They analyzed these lines in order to check the accuracy of the 
oscillator strengths, by comparing the abundance obtained with these lines and the OH 
lines of the A v = 1 sequence. Bonnell & Bell (1993) used the HM74 model atmosphere and 
obtained an abundance of = 8.76 (cr = 0.09 dex) with the first overtone vibration-rotation 
OH hues, that is 0.12 lower than the value they found with the fundamental OH hues. After 
considering the mean difference in ^/-values and equivalent widths with respect to our work, 
their result is Aq — 8.71 ± 0.09 dex, which is lower than the abundance we found with the 
HM74 model atmosphere {Aq = 8.80 dex ± 0.06 dex). M. Asplund have also calculated the 
abundance of the first-overtone OH lines using the HM74 model atmosphere and a subset 
of five lines, resulting in Aq = 8.79 ± 0.06 dex (after normalizing his result with the 15 OH 
lines). Considering systematic errors of about 0.02 dex, and mainly considering the large 
statistical errors involved, the results with the HM74 model atmosphere are in reasonable 
agreement. 

Before comparing our results with those obtained with other lines analyzed with 3D 
hydrodynamical atmospheres, we estimate the total error of our oxygen abundance for the 
3D case. The quadratic sum of the observational error (fa 0.06 dex), the error of the gf- 
values (Ri 0.03 dex), and a systematic error of (at least) 0.02 dex, results in a total error of 
5Ao = 0.07 dex. 

4. Comparison with 3D studies 

Considering that Kurucz model atmospheres are widely used in the literature, we have 
also employed a Kurucz solar model to determine oxygen abundances using the [O 1], O 1 
and Av = 0, 1 OH lines. For comparison purposes, the abundance from these spectral 
features was also obtained using the <3D> solar model. In short, the equivalent widths and 
^/-values of the [O I] and the O 1 777 nm triplet was obtained from A04 (note that Wx for 
the forbidden lines are cleaned from the Ni 1 and CN blends), and the equivalent widths for 
the Av — 1 and pure rotation OH lines were measured from the WL03 and ATMOS solar 
atlas, respectively. The gf -values for the molecular lines are from G98, which are essentially 
identical to the /-values employed by A04. NLTE corrections for the O I triplet were taken 



- 7- 



from A04. 

In Table 2 a comparison of our results with those obtained by A04 is shown. For the 
3D case we can see that the molecular lines seem to give shghtly lower abundances than the 
[O I] and O I lines. This has been already noted by AP04, who suggested that this difference 
could signal that the results derived from the infrared OH lines are less rehable than those 
obtained from the forbidden and permitted lines. 

Are these small differences in the abundance obtained from different features due to 
errors in the transition probabilities? The gi/-values of the pure rotation OH lines seem 
highly reliable. In fact, an accurate dipole moment (1.6676 ± 0.0009 Debye) for v — 0, 
J = 9/2, have been measured by Meerts &: Dynamus (1973). As the vibration-rotation 
interaction of the pure rotational lines is very small, an extrapolation of the DMF (using 
theoretical calculations) for lines with higher J and v seems safe (see e.g. Fig. 2 of Sauval et 
al. 1984). Goldman et al. (1983) adopted a constant DMF for their analysis of the oxygen 
abundance in the Sun from the pure rotation OH lines. They made a small correction to the 
experimental value derived by Meerts & Dynamus (1973), to take into account the higher J 
values of the solar lines. In Fig. 3a we show the difference between the gf -values derived by 
G98 and Goldman et al. (1983), as can be seen, neglecting the variation of the DMF results 
in an error of less than 1%. A comparison between the (^/-values obtained from G98 and 
Goorvitch et. al. (1992) is also shown is Fig. 3a. Goorvitch et al. adopted Nelson's DMF, so 
the agreement with G98 is not surprising. The error bars shown in Fig. 3a. were estimated 
from the maximum error that is made if a constant DMF were adopted (for typical solar 
pure rotation OH lines with v = 0, 1, 2, 3). An error of 0.01 dex for the 5f/-values of the Av 
— sequence is not unreasonable. 

The errors for the fundamental {Av= 1) OH lines were estimated by N90, as discussed 
in §3. In Fig. 3b are shown these errors and the difference between the gf -values derived 
from G98 and N90. In a different approach, G98 made error estimates by making slight 
modifications to the shape of the experimental DMF by N90. Using these results we obtained 
the dashed error bars shown in Fig. 3b. The typical error for the ^'/-values of the Av — 1 
hues seem to be about 0.035 dex. 

In §3 we discussed the errors for the transition probabihties of the Av — 2 OH lines, 
they are shown in Fig. 3c, and are about 0.03 dex. 

APOl states that the uncertainty of the (//-values of the forbidden lines is about 0.02 
dex. Lambert (1978) suggested log gf = -9.75 for the 6300 A line, 0.03 dex lower than the 
gf value adopted by APOl and A04, and the NIST database of critically selected transition 
probabilities suggests log gf — -9.78, 0.03 lower than the value recommended by Lambert 
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(1978). We adopt an uncertainty of 0.03 dex in the /-value of the forbidden hues. Another 
important uncertainty in the O abundance derived from the 6300 and 6363 A [O I] hues is due 
to blends with Ni I and CN lines. In order to evaluate these uncertanties, we have computed 
synthetic spectra of the CN lines that are blended with the [O I] lines. The contribution of 
CN lines is negligible for the 6300.3 A line, but it is very important for the 6363.7 A line. As 
already noted by A04, the main blend is due to the (10,5) Q2 25.5 CN line. Our calculation 
suggests a similar correction to that proposed by A04 (0.5 mA). Fortunately, it is possible to 
observationally check this result, by measuring the (10,5) Ql 25.5 CN line, which obviously 
has the same excitation potential and gf-value as the Q2 CN line. The Ql CN hue at 6365.0 
A is apparently unblended, but the line is very weak and the continuum is difficult to define. 
After trying several solar atlases, we found that this line was better defined in the Wallace 
et al. (1998) FTS solar spectrum at disk center. When a correction to solar fiux is made^, 
Wx = 0.35 mA is obtained^, which is lower than the 0.5 mA employed by A04. The error in 
Ao[6363] due to this uncertainty in the CN blend is 0.04 dex. On the other hand, adopting 
a 20% uncertainty for VI/;^(Ni I), results in an error of about 0.03 dex in Aq determined from 
the 6300 A line. These errors from the CN and Ni I blends have been considered in column 
5 of Table 3. 

For the permitted lines, A04 used transition probabilities from the NIST database. The 
difference between these oscillator strengths and two different calculations given by Biemont 
et al. (1991) is 0.025 dex, so we adopt 0.03 dex as the error in the gf -values of the permitted 
fines. 

The adopted errors are shown in Table 3, where are also shown the statistical uncertain- 
ties (cr), the errors due to NLTE effects, and estimates of other errors made in the modeling 
(see Table 3), as well as the quadratic sum of the errors. The error due to NLTE effects is 
only important for the O I lines. In particular, for the 777 nm triplet, the NLTE correction is 
A ^ —0.21 dex, but it depends on the empirical correction factor Sh- The NLTE correction 
could be 0.03 dex lower or higher, if Sh = (no collisions with H) or Sh = 1 (Drawin-like 
formula) is adopted, respectively (AP04). Conservatively, we have adopted an NLTE error 
of 0.04 dex, slightly higher than this 0.03 dex uncertainty, and for the abundance obtained 
from the permitted lines we will consider the mean abundance obtained by A04 and AP04. 

In Fig. 5 we show the oxygen abundances obtained from the different lines in the 3D 
analysis (Table 2), and the error bars obtained above (Table 3). The recommended solar 



^A04 used the Kurucz et al. (1984) solar flux atlas, instead of a disk center spectrum 

^Note that if the true continuum is significantly higher than the local continuum, then this measurement 
is a lower limit 
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oxygen abundance by A04 {Aq — 8.66 ± 0.05) is also shown. As can be seen, within the 
errors, all different spectral features, from the optical to the far infrared, are in good agree- 
ment. The weighted- average (using the inverse square of the errors as weight) considering 
all the lines is only 0.02 dex lower than the value recommended by A04. 

5. Summary and Conclusions 

A low oxygen abundance of Aq ~ 8.6 dex (3D, <3D> and MARCS models) was obtained 

from extremely weak (log {Wx/Xj ~ —7.2) first-overtone vibration-rotation OH lines present 
in the H band. These lines are about a factor of 5 weaker than the fundamental vibration- 
rotation OH lines present in the L band of the solar spectrum. The OH lines used in this 
work arc the weakest ones ever used to obtain the oxygen abundance of the Sun (e.g. for log 
{Wx/X) = -7.2, an optical line at 6000 A has Wx just 0.4 mA). 

The present results support the low solar oxygen abundance recently obtained from other 
spectral features using 3D hydrodynamical model atmospheres (as summarized in Table 2). 

Comparing the results obtained with different spectral features and model atmospheres, 
we can see that the best agreement with the 3D simulations is obtained with the <3D> 
model, this is not unexpected, since this model represents the temporal and spatial average 
of different 3D calculations. As stressed by A04, the <3D> model cannot reproduce all 3D 
effects, due to the lack of atmospheric inhomogeneities; this can be seen comparing columns 
<3D> and 3D of Table 2. Note that the first overtone OH hues are the less sensitive of 
the OH hues to the 3D effects, probably because these hues are formed deeper than the 
fundamental and pure rotation OH lines. 

This new determination of the oxygen abundance in the Sun is consistent with the 
abundance obtained in B and late O stars (Daflon et al. 2003), solving the problem of the 
Sun being more oxygen-rich than the ISM. Nevertheless, the depth of the solar convection 
zone computed with the new solar abundance and OPAL opacities (Bahcall & Pinsonneault 
2004) is in conflict with the measured value of the base of the convection zone obtained 
by hehoseismology (Basu & Antia 1997, 2004). Accordingly to Basu & Antia (2004), the 
problem could be solved with an upward revision of either the solar abundances or the opacity 
tables. Bahcall et al. (2004a,b) have estimated that an increase of about 10% in the OPAL 
opacity tables is required to solve the "convective zone problem" and to bring the computed 
sound speed profile into agreement with helioseismological observations. Recently, Seaton & 
Badnell (2004) have shown that new calculations from the Opacity Project (OP) result in 
opacities about 5% larger than those given by OPAL (for temperatures and densities near 
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the base of the solar convective zone), alleviating the difference between solar models and 
helioseismology, and bringing support for a low solar oxygen abundance. 

I thank the anonymous referee for his valuable suggestions and noting recent pre-prints 
on helioseismology and OP opacities; M. Asplund for his comments, sending the MARCS 
solar model and performing 3D calculations; L. Wallace for providing a printed copy of the 
NSO solar atlas; and J. G. Cohen for her comments. NSO/Kitt Peak FTS data used here were 
produced by NSF/NOAO. I am grateful for partial support from NSF grant AST-0205951 
to J. G. Cohen. 
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Fig. 1. — Observed solar spectrum (WL03) around 1.69 //m. The vertical lines indicate three 
OH hues used in this work. The strongest atomic features are identified. The profile of the 
OH line at 16904.3 A is shown in the inset, where only the upper 0.7 % of the spectrum is 
plotted. 




Fig. 2. — Comparison between our Wx measurements and those by Bonnell & Bell (1993). 
The hne depicts perfect agreement. 
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Fig. 3. — Differences between the (^/-values of infrared OH lines calculated by G98 and (a) 
Goorvitch et al. (1992) (filled circles) and Goldman et al. (1983) (open circles); (b) and 
(c) N90 (filled circles). The error bars given in (a) were obtained supposing that the DMF 
is constant. The sohd error bars in (b) and (c) were estimated by N90 taken into account 
the error of its experimental DMF. The dashed error bars in (b) were estimated from slight 
variations to the Nelson's DMF, as given in the last column of Table 3 by G98. 
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Fig. 4. — Solar oxygen abundances obtained with the <3D> model and first overtone infrared 
OH lines (filled circles) vs. (a) excitation potential, (b) W\/\, and (c) wavelength. The mean 
abundance (Aq = 8.59 dex) is shown by dotted lines, and the a (0.06 dex) error bar is shown 
in the middle panel. 
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Fig. 5. — Solar oxygen abundance derived from different lines in 3D analyses (Table 2). The 
error bars are from Table 3. The recommended value given by A04 [Aq = 8.66 ± 0.05 dex) 
is shown by the dotted and dashed hnes. 
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Table 1. Equivalent widths of OH lines 



(A) 






ID^ 




Xexc 

(eV) 


(d 


5 gf 
ex) 


Wx 
(ml) 


15278 


52 


(2,0) 


Pic 


9.5 


0.205 


-5 


382 


0.9 


15409 


17 


(2,0) 


P2c 


9.5 


0.255 


-5 


365 


1.1 


15568 


78 


(2,0) 


Pie 


11.5 


0.299 


-5 


269 


1.2 


16052 


76 


(3,1) 


Pie 


9.5 


0.639 


-4 


910 


1.25 


16192 


13 


(3,1) 


P2c 


9.5 


0.688 


-4 


893 


1.4 


16368 


13 


(3,1) 


Plf 


11.5 


0.731 


-4 


797 


1.3 


16456 


04 


(2,0) 


Plf 


16.5 


0.609 


-5 


048 


0.9 


16534 


58 


(3,1) 


Pie 


12.5 


0.782 


-4 


746 


1.35 


16605 


46^= 


(4,2) 


Pic 


7.5 


0.982 


-4 


755 


0.75 


16655 


99 


(2,0) 


Pie 


17.5 


0.682 


-5 


010 


1.15 


16872 


28^ 


(2,0) 


Pie 


18.5 


0.759 


-4 


975 


0.7 


16886 


28 


(4,2) 


P2c 


8.5 


1.059 


-4 


662 


0.8 


16904 


28 


(3,1) 


Pic 


14.5 


0.897 


-4 


654 


1.0 


16909 


29 


(3,1) 


Plf 


14.5 


0.898 


-4 


654 


1.3 


17618 


89 


(4,2) 


Pie 


13.5 


1.244 


-4 


403 


1.3 



"^identifications in the solar spectrum by WL03 

''(v',v") branch J" 

'^lines identified in this work 
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Table 2. Solar oxygen abundance from different spectral features and model atmospheres. 



lines 3D HM74 MARCS <3D> Kurucz 



[O I]^ 
OP 

OH (Aw = 0)^ 
OH (A^; = l)'^ 
OH {Av = 2)^ 



8.68 ± 0.01 
8.67= ± 0.02 
8.65 ± 0.02 
8.61 ± 0.03 
8.57 ± 0.06 



8.76 ± 0.02 

8.64 ± 0.08 
8.82 ± 0.01 
8.87 ± 0.03 
8.80 ± 0.06 



8.72 ± 0.01 

8.72 ± 0.03 
8.83 ± 0.03 
8.74 ± 0.03 
8.61 ± 0.06 



8.75 ± 0.01 

8.68 ± 0.03 
8.74^ ± 0.05 

8.69 ± 0.06 
8.59 ± 0.06 



8.79 ± 0.01 

8.67 ± 0.01 
8.80'= ± 0.06 
8.82 ± 0.05 
8.71 ± 0.06 



Note. — The quoted uncertainties are only the a (line-to-line scatter) 

Results from A04 for the 3D, HM74 and MARCS models, and abundances obtained in 
this work for the <3D> and Kurucz models 

''This work, except for the 3D case. The 3D calculations were performed by M. Asplund, 
who obtained 8.58 ± 0.05 for 5 OH lines of the Av = 2 sequence, using this result we found 
Ao = 8.57 ± 0.06 for the 15 OH lines 

■^This is the mean of the abundance obtained by A04 and AP04 

'^Ao = 8.75 ± 0.03 dex is obtained when only lines with Wx = 15-100 mA are employed 
''Ao = 8.82 ± 0.03 dex for the lines with Wx = 15-100 mA 



Table 3. Errors in the oxygen abundance (dex) 



lines 




a'' 


gfs 


NLTE 


others'' 


total 


[01] 




0.01 


0.03 


0.00 


0.04 


0.05 


01 




0.02 


0.03 


0.04 


0.02 


0.06 


OH (At; = 


0) 


0.01 


0.01 


0.00 


0.05 


0.05 


OH {Av = 


1) 


0.03 


0.04 


0.00 


0.02 


0.05 


OH {Av = 


2) 


0.06 


0.03 


0.00 


0.02 


0.07 



^line-to-line scatter 

''minimum uncertainties in the 3D modeling, due to 
systematic errors in the continuum opacities, partition 
functions, etc. For the [O I] lines errors due to blends 
with Ni I and CN lines arc considered. A higher error 
for the pure rotation lines was adopted because the 3D 
modeling is less successful for these lines (see Fig. 10 of 
A04) 



